>> x=[1985:1:2003]';
>> X=[ones(19,1) x];
>> Y=[470 481 492 513 571 597 630 694 708 719 720 750 770 801 795 828 862 937 1006]';
>> [b,bint,r,rint,stats]=regress(Y,X);
>> b,bint,stats

b =

  1.0e+004 *

   -5.3283
    0.0027


bint =

  1.0e+004 *

   -5.8204   -4.8361
    0.0025    0.0030


stats =

    0.9692  535.6016    0.0000  780.0594

>> rcoplot(r,rint)
>> z=b(1)+b(2)*x

z =

  458.6526
  485.7263
  512.8000
  539.8737
  566.9474
  594.0211
  621.0947
  648.1684
  675.2421
  702.3158
  729.3895
  756.4632
  783.5368
  810.6105
  837.6842
  864.7579
  891.8316
  918.9053
  945.9789

>> plot(x,Y,'k +',x,z,'r')
>> rcoplot(r,rint)
>> 







>> x=[64 70.3 76.9 81.5 90.3 153 180.9 200.8 215.5 218.3 249.4 251.6];
>> y=[0.119 0.1219 0.4321 0.1200 0.0663 0.0185 0.0351 0.0224 0.0107 0.0196 0.0068 0.0254];
>> A=polyfit(x,y,2)

A =

    0.0000   -0.0033    0.3791

>> z=polyval(A,x);
>> plot(x,y,'k +',x,z,'r')
>> 